Higher Twists and Color Polarizabilities in the Neutron 
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O ■ Abstract 

-xt- : 

Color polarizabilities of the neutron are extracted from data on the lowest moment of the spin- 

^t- ; 

dependent g\ structure function. New data in the resonance region from Jefferson Lab at Q 2 < 

r-j ■ 

1 GeV 2 , in combination with world data at higher Q 2 , allow a systematic determination of the 

Ph. 1/Q 2 corrections, and provide the first constraints on l/Q 4 corrections. The results suggest that 

r- 1 \ 

higher-twist effects in the neutron are small, and that quark-hadron duality may be approximately 

> ■ 

valid, even down to Q 2 ~ 1 GeV 2 . 
PACS numbers: 12.38.Aw, 12.38.Qk, 13.60.Hb 
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The nature of quark confinement in QCD remains a central mystery of strong interaction 
physics. Any satisfactory solution to this problem requires understanding of the dynamics of 
quark and gluon interactions at large distances, where QCD becomes highly nonperturbative. 

The properties of the quark and gluon constituents of the nucleon are usually determined 
by examining how the quark and gluon fields respond to external probes. Deep inelastic elec- 
tron scattering at high energy transfer v and momentum transfer squared Q 2 has provided a 
wealth of information on the single-particle momentum and spin distributions of quarks and 
gluons in the nucleon. Interactions between quarks and gluons at short distances, which give 
rise to scaling violations, can be described perturbatively through the Q 2 evolution equa- 
tions. To probe directly long-range quark-gluon interactions, on the other hand, requires 
the study of structure functions at intermediate values of Q 2 , in the pre-asymptotic region. 

In this region a transition takes place between deep inelastic scattering at high Q 2 , where 
structure functions are determined by incoherent scattering from partons, and the low Q 2 
region which is characterized by resonance dominance of structure functions [jj. The re- 
gion of intermediate Q 2 (Q 2 ~ 1 GeV 2 ) exposes quark-gluon correlation effects, which are 
otherwise suppressed in the scaling regime, while still allowing an expansion in 1/Q 2 to 
be meaningful. An example of such effects is the so-called color polarizabilities (not to be 
confused with the more familiar electromagnetic polarizabilities), which describe how the 
color electric and magnetic fields respond to the spin of the nucleon. 

In this paper we analyze data from Jefferson Lab Hall A experiment E94-010 2] on the 
neutron spin dependent structure function g\ at Q 2 < 1 GeV 2 , which when consistently 
combined with the world high Q 2 data, enable a first accurate determination of the color 
polarizabilities in the neutron. In addition, we explore the closely related phenomenon of 
quark-hadron duality, where the gi spin structure function in the resonance region averages 
to the scaling function measured at high Q 2 . This dualityhas been explored in recent 
experiments for the proton and deuteron structure functions p|, |^ , but has not been studied 
in the case of neutrons. 

The analysis makes use of the operator product expansion (OPE) in QCD, where at large 
Q 2 the lowest moment, T"(Q 2 ), of the spin structure function gi of the neutron is expanded 
in inverse powers of Q 2 , 




(1) 



with the coefficients /x™ related to nucleoli matrix elements of operators of twist < r. Here 
twist is defined as the mass dimension minus the spin of an operator, and x = Q 2 /2Mv is 
the Bjorken x variable, with M the neutron mass. For each twist the Q 2 dependence in /i™ 
can be calculated perturbatively as a series in a s . Note that the application of the OPE 
requires summation over all hadronic final states, including the elastic at x = 1. 

The leading-twist (twist-2) component, is determined by matrix elements of the axial 
vector operator ^7^75^, summed over various quark flavors. It can be decomposed into 
flavor triplet («7a), octet (a§) and singlet (AS) axial charges, 

/4(Q 2 ) = C ns (Q 2 ) (~9a + ^a 8 ) +C s (g 2 )^AS , (2) 

where C ns and C s are the nonsinglet and singlet Wilson coefficients The nonsinglet 
triplet axial charge is obtained from neutron /3-decay, qa = 1.2670(35) p, while the octet 
axial charge is extracted from hyperon weak decay matrix elements assuming SU(3) flavor 
symmetry, a§ = 0.579(25) [6J. In order to factorize all of the Q 2 dependence into the Wilson 
coefficients, we use the renormalization group invariant definition of the matrix element of 
the singlet axial current, AS = AS(Q 2 = 00). 

The higher-twist contribution to T™(Q 2 ) can be obtained by subtracting the leading- twist 
term from the total, 

ATUQ 2 ) = TUQ 2 )-^(Q 2 ) = ^f 1 + ^f 1 + o (±) . (3) 

The coefficient of the 1/Q 2 term contains a twist-2 contribution, a 2 , and a twist-3 term, d 2 , 
in addition to the genuine twist-4 component, f 2 

ft = l -M 2 (a n 2 + Ad n 2 + 4#) . (4) 

The twist-2 term a 2 arises from the target mass correction Q], and is related to the second 
moment of the twist-2 part of the g™ structure function. 

Interaction terms appear through the twist-3 and twist-4 contributions. The coefficient 
d 2 is given by the twist-3 part of the x 2 -weighted moment of the g 2 structure function, and 
corresponds to a matrix element of an operator which involves both quark and gluon fields 

0. 

The twist-4 contribution to /i^ is defined by the matrix element 

M 2 S^ = e " ^9 & ^\ N ) > ( 5 ) 



where G^ u is the dual gluon field strength tensor, and g is the strong coupling constant. 

The twist-3 and 4 operators describe the response of the collective color electric and 
magnetic fields to the spin of the neutron. Expressing these matrix elements in terms of the 
components of G^ v in the nucleon rest frame, one can relate d 2 and f£ to color electric and 
magnetic polarizabilities. These are defined as 
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X n E 2M 2 S = (N\ ja x E a \N) , (6) 

X n B 2M 2 S = (N\ f a B a \N) , (7) 

—* — — * — * 

where S is the neutron spin vector, j£ = —gil)^H a ip is the quark current, and E a and B a 

are the color electric and magnetic fields, respectively. In terms of d 2 and f 2 the color 
polarizabilities can be expressed as 

X n E = \ (2d n 2 + ) , x n B = \ M - ft) . (8) 



An analysis of the twist-4 matrix elements was carried out in Ref. |rjj using data on 
from the E143 experiment at SLAC Il2ll at Q 2 = 0.5 and 1.2 GeV 2 . Subsequently 

13] and E155 [14], the HERMES Collaboration at DESY |lj, 
16| | have collected data on T™ over a larger range of Q 2 . The 



SLAC experiments E154 
and the SMC at CERN 
precision of the earlier data in the resonance region was rather poor, and no data existed for 
0.5 < Q 2 < 1 GeV 2 , which precluded the size of the twist-4 matrix elements from being well 
constrained. More recently, precision data from Jefferson Lab Hall A experiment E94-010 
Q] were collected on at several Q 2 values below 1 GeV 2 . These data allow one for the 
first time to determine systematically the size of the f£ matrix element, and estimate the 
l/Q 4 correction. 

Before proceeding with the extraction, it is important to ensure that the various data sets 
are consistent with each other, in terms of the assumptions and extrapolations beyond the 
measured region employed in each analysis. For instance, each of the experimental analyses 
uses a different extrapolation to x — 0, which makes a naive comparison of the published 
T™ values problematic. To this end we have reexamined all of the previous data on T™ using 
the same set of assumptions for the x — > as well as the x — > 1 extrapolations. 

To evaluate for each experiment, we added to the quoted r" value from the mea- 
sured region a contribution of the unmeasured \ow-x region using the Bianchi-Thomas (BT) 



parametrization 



17| up to W 2 = 1000 GeV 2 . This method was applied to the SLAC 



experimental results jl^l . Q]. In the case of the HERMES and Jefferson Lab exper- 
iments, the published T™ values already include an estimate from the unmeasured low- 
x region using the BT parametrization, while for the SMC experiment data exist above 
W 2 = 1000 GeV 2 . Except for the HERMES data, the remaining contribution down to 
x = was then evaluated assuming g\ to be constant: gi(x,Q 2 ) = gf T (x min , Q 2 ) for the 
SLAC and JLab experiments, with x min defined for each Q 2 by W 2 = 1000 GeV 2 , and 
<(x,Q 2 ) = gf MC (x^ = 0.003, Q 2 = 10 GeV 2 ) for the SMC experiment. Since the SLAC 
E155 Collaboration |lj] only presents a moment deduced from a next-to-leading order (NLO) 
analysis using world data, we do not include this result in our fit. Nevertheless, the NLO 
result for at Q 2 = 5 GeV 2 including the SLAC E155 data is very close to the result of 
the SLAC E154 NLO analysis Q- 

To include the unmeasured x — > 1 region, we used a new fit of the world neutron data 
which includes the high-precision A, data of the M OT Lao E99-U7 e„nt Q at 
large x. Since the elastic contribution is included separately, the maximum value of x is 
defined for each experiment by the pion electroproduction threshold. The resulting total 
moments from the world data are plotted in Fig. 1 for 0.5 < Q 2 < 10 GeV 2 , where 
the total uncertainty in each data set is the quadratic sum of the statistical and systematic 
uncertainties. The Jefferson Lab experiment E94-010 (filled circles) extends the range of Q 2 
with precision data below Q 2 = 1 GeV 2 . In all cases the data include both the inelastic and 
elastic contributions, with the latter taken from the fit in Ref . j^j] . 

The twist-2 contribution fi^ is determined by fitting the neutron data in Fig. 1 assuming 
there are no higher twists beyond Q 2 = 5 GeV 2 , from which we obtain AS = 0.35 ± 
0.08, where the uncertainty is statistical. Using this central value, the twist-2 contribution 
is illustrated in Fig. 1 by the shaded band, with the extrema representing the range of 
uncertainty associated with the value of a s in the Wilson coefficients. The exact value of 
AS depends somewhat on the x — > behavior assumed in the extrapolation beyond the 
measured region. However, since the higher-twist contributions are determined from the 
relative variation in T™ from high to low Q 2 , the absolute normalization of the leading- twist 
contribution does not play a major role in determining f£. 

The higher-twist contribution AT", obtained by subtracting the leading- twist curves in 
Fig. 1 from data on the total moment T", is shown in Fig. 2 as a function of 1/Q 2 for 
AS = 0.35. Here we have used = —0.0031(20) for the target mass corrections, obtained 
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FIG. 1: Q 2 dependence of T™ from various experiments. The error bars are a quadratic sum of 
statistical and systematic uncertainties. The twist-2 contribution from Eq. @ is given by the 
band with AS = 0.35, and its width represents the uncertainty in a s . The elastic contribution is 
indicated by the dashed curve. 



from a fit to the world neutron data |20j at Q 2 = 5 GeV 2 , and the value d? = 0.0079(48) for 
the twist-3 matrix e, m ent obtained f r0m SLAC expert E15 5X Q. At this Q 2 va.ue 
and c?2 are dominated by their leading-twist contributions. 

While the Q 2 evolution of the (twist-2) is straightforward, the evolution of higher-twist 
structure functions is in general rather more involved. For the twist-4 matrix element 
the Q > evoMion w a8 counted in Refs. Q Q to leading logar ithnne oj, In this an^sis 
we assume the leading-twist values for and d 1 ^ at Q 2 = 5 GeV 2 and use the results from 
Refs. [3, 13 to account for the logarithmic Q 2 dependence of f£. In practice, the inclusion 
of a s dependence of the 1/Q 2 corrections has very little influence on the values of the higher 
twists that we extract. 

The solid curve in Fig. 2 represents a 2-parameter minimum \ 2 fit to the Ar™ data for 
Q 2 > 0.5 GeV 2 , using Eq. (jHJ) with and the l/Q 4 correction /Xg as free parameters. We 
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FIG. 2: Higher-twist correction Ar™ versus 1/Q 2 . The world data points include statistical (inner 
ticks) and total uncertainties (outer ticks), except for those of HERMES and JLab E94010, for 
which only statistical uncertainties are shown with error bars, with systematic uncertainties indi- 
cated by the dark bands at the bottom of the figure. The solid curve is a 2-parameter (fg and 
//g) fit to the Q 2 > 0.5 GeV 2 data, while the dashed curve is a 1-parameter (fg only) fit to the 
Q 2 > 1 GeV 2 data. The band around the solid curve represents the uncertainty of the fit due to 
statistical uncertainties, and the light band at the bottom of the figure corresponds to the total 
uncertainty. 

neglect any possible Q 2 dependence in /i 6 itself, which should be a reasonable assumption 
within the present uncertainties. The best fit values for the twist-4 and 1/Q 4 corrections, 
using only the statistical uncertainty for each experiment, are found to be 

/ 2 n = 0.033 ± 0.005 , = (-0.019 ±0.002)M 4 , (9) 

normalized at Q 2 = 1 GeV 2 . Including the total systematic uncertainty for each experiment, 
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we find 

/ 2 n = 0.034 ± 0.043 , $ = (-0.019 ±0.017)M 4 . (10) 

With this latter value of fg, the overall l/Q 2 correction to the leading-twist term, at a scale 
Q 2 = 1 GeV 2 , including the target mass, twist-3 and twist-4 contributions extracted from 
data, is fi^ = (0.019 ± 0.024)M 2 . Combining this with the l/Q 4 term, the total higher-twist 
contribution Ar™ is almost exactly zero at Q 2 = 1 GeV 2 . 

The reliability of truncating the expansion at l/Q 4 order down to Q 2 = 0.5 GeV 2 can 
be tested by adding a /ig/Q 6 term. Performing a 3-parameter fit results in almost identical 
values for fg and /ig, with /ig = (0.00 ± 0.03)M 6 . On the other hand, if the data are 
fitted with just the l/Q 2 term, one generally finds a smaller value for fg, namely f£ = 
—0.014 ± 0.010, with a somewhat larger x 2 ( see a ls° Ref- 24|). This clearly illustrates the 
influence of the JLab E94-010 data, which suggest a negative trend in Ar™ for Q 2 < 1 GeV 2 . 
Excluding the JLab results and fitting only the Q 2 > 1 GeV 2 data, the uncertainty on fg 
becomes 2-3 times larger, and the ^ coefficient essentially unconstrained. A 1-parameter 
fit to the Q 2 > 1 GeV 2 data, indicated by the dashed curve in Fig. 2, gives a slightly smaller 
/ 2 n value, / 2 n = 0.012 ± 0.029. 

The extracted value of f£ can be compared with the earlier analysis in Ref. [llj . which 
found f£ = 0.07 ± 0.08 using a 1-parameter fit to Q 2 > 1 GeV 2 data. The present analysis 
improves the uncertainty on f% by a factor ~ 3 using the non-JLab data, and by a factor ~ 8 
including the JLab data. Of course, the uncertainty on / 2 in Ref. [llj would be significantly 
larger had a 2-parameter fit been employed. 

The value of f£ extracted from the data can be compared with nonperturbative models. 
For example, QCD sum rules Js[ |^ yield negative values for the twist-4 matrix elements: 



^Isummic = -0.018(17) and -0.013(6) from Refs. [25jand jj, respectively. On the 
other hand, calculations within the instanton model give [26j ,/2l instanton = 0.038, in good 
agreement with our central value. Estimates of the r = 4 matrix elements using the MIT 
bag model [ll]], evolved from the bag scale, give /^Ibag = 0-0, which is also consistent with 
our result. 

Combining the extracted fg values with the value obtained by the E155X Collaboration 



22|, the color electric and magnetic polarizabilities, normalized at Q 2 = 1 GeV 2 , are found 
to be 

X n E = 0.033 ± 0.029 , Xb = -0.001 ±0.016 . (11) 



These results indicate that both the color electric and magnetic polarizabilities in the neutron 
are relatively small, with the magnetic polarizability consistent with zero. The small values 
of the higher-twist corrections suggest that the long-range, nonperturbative interactions 
between quarks and gluons in the neutron are not as dominant at values of Q 2 > 0.5 GeV 2 as 
one may have expected. This would imply strong cancellations between neutron resonances 
resulting in the dominance of the leading-twist contribution to T™. The results presented 
here therefore provide a spectacular confirmation of quark-hadron duality in spin- dependent 
structure functions. A dedicated experiment [^J to study duality in neutron spin structure 
functions in the resonance region is currently being analyzed. 

In summary, we have determined the first moment of the neutron spin structure function 
gi from the world data using consistently the same extrapolation method to x = and 
x = 1. From these data we extracted the size of the twist-4 matrix element f£ of the neutron. 
The neutron color polarizabilities were determined by combining f£ and the twist-3 matrix 
element from SLAC experiment E155X, which leads to a small and slightly positive value 



the g\ and g% structure functions at Q 2 ~ 1 — 4 GeV 2 will reduce the uncertainty in the 
extracted higher-twist coefficients, as will better knowledge of a s in the intermediate Q 2 
region. 
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